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Equation of state for chromium nitride has been debated in the literature in connection with a 
proposed collapse of its bulk modulus following the pressure induced transition from the paramag- 
netic cubic phase to the antiferromagnetic orthorhombic phase [F. RivaduUa et ai, Nat Mater 8, 974 
(2009); B. Ailing et al, Nat Mater 9, 283 (2010)]. Experimentally the measurements are complicated 
due to the low transition pressure, while theoretically the simulation of magnetic disorder represent 
a major challenge. Here a first-principles method is suggested for the calculation of thermody- 
namic properties of magnetic materials in their high temperature paramagnetic phase. It is based 
on ab-initio molecular dynamics and simultaneous redistributions of the disordered but finite local 
magnetic moments. We apply this disordered local moments molecular dynamics method to the 
case of CrN and simulate its equation of state. In particular the debated bulk modulus is calculated 
in the paramagnetic cubic phase and is shown to be very similar to that of the antiferromagnetic 
orthorhombic CrN phase for all considered temperatures. 



I. INTRODUCTION 



Chromium nitride is a material which combines practi- 
cal and industrial relevance as a component in protective 
coatings''" with fascinating fundamental physical phe- 
nomena. The latter include a phase transition with a 
magnetically driven lattice distortion' between an anti- 
ferromagnetic orthorhombic low temperature phase and 
a paramagnetic cubic high temperature phase . The im- 
portance of strong electron correlations as well as the 
necessity to model the paramagnetic state using finite 
disordered local moments have been recently shown''*'. 
Important issues, such as the impact of the phase tran- 
sition on the compressibility of the material ' as well as 
on the electrical conductivity'^''" are still subjects of an 
intense discussion. 



The core problem of obtaining a complete understand- 
ing of these phenomena and properties on the most fun- 
damental level of physics arises from the difhculty of sim- 
ulating the paramagnetic high-temperature phase from 
first principles. In this work we first discuss the method- 
ologies that have been used in theoretical treatments of 
paramagnetism. Then we present a practical scheme for 
calculating thermodynamic properties, in particular the 
equation of state, of a paramagnetic material at elevated 
temperature merging ab initio molecular dynamics (MD) 
and the disordered local moments model (DLM). This 
DLM-MD technique is then applied to investigate the 
infiuence of temperature and pressure on the compress- 
ibility of CrN. We show that the change of the bulk mod- 
ulus of CrN upon the pressure induced phase transition 
is minimal, strengthening conclusions from earlier static 
calculations'''^ which questioned its reported collapse' . 



II. MODELING THE PARAMAGNETIC STATE 
A. Background 

A theory that describes the finite-temperature aspects 
of itinerant electron magnets has to take into account 
the existence of local magnetic moments present above 
the magnetic transition temperature, the Curie tempera- 
ture Tc or the Neel temperature for a ferromagnetic 
or an antiferromagnetic material respectively . At the 
same time, the majority of methods used for ab initio 
electronic structure calculations nowadays are based on 
the density functional theory (DFT) in the local (local 
spin density, LSDA) or semi- local (generalized gradient, 
GGA) approximations. While they are known to give 
an accurate description of the ground state properties 
of magnetic systems^^", its straightforward generalization 
to finite temperatures leads to quantitative, as well as 
qualitative errors' Indeed, Tc of transition metals are 
overestimated by a factor of five and there are no mo- 
ments and no Curie- Weiss law above Tc- A solution to 
this problem should in principle be sought in the physics 
of strongly correlated electron systems. In particular, 
the dynamical mean-field theory (DMFT) , combined 
with LDA band structure calculations has been used with 
success for simulations of finite-temperature magnetism 
in Fe and Ni^"'. However, its application to the study 
of structural phase transition in Fe^'' had to neglect a 
contribution from the lattice dynamics, because of pro- 
hibitively high computational cost and difficulties in cal- 
culating forces between atoms . 

At the same time, it is realized that LSDA calculations 
at zero temperature can provide valuable information for 
the description of the finite temperature magnetism. One 
way of doing this is to extract magnetic interactions in 
the form of exchange constants for a classical Heiscn- 
bcrg model or magnetic "forces" (the first variation of 
the total energy for a differential rotation of a local mo- 
ment)^"'"" from DFT calculations and to use them in 
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statistical mechanics^^"^^ or in spin-dynamics'^*^''"" simu- 
lations of magnetic properties at elevated temperatures. 
Another useful approach is given by the so-called Disor- 
dered Local Moment model, introduced by Hubbard 
and Hasegawa ' ' ' ■ , and combined with the LSDA-DFT 
by Gyorffy et al.^ ''. Within the DLM picture, the lo- 
cal magnetic moments exist in the paramagnetic state 
above the magnetic transition temperature, but are fully 
disordered. The magnetically disordered state can be de- 
scribed as a pseudo-alloy of equal amounts of atoms with 
spin up and spin down orientations of their magnetic mo- 
ments, and its electronic structure and the total energy 
can be calculated within the conventional alloy theory 
using the coherent potential approximation (CPA) or 
the supercell technique''. The methodology is highly suc- 
cessful in applications to many materials systems, rang- 
ing from steels'^ '' ' ^ to actinides * ', and its generalization 
for the case of partial magnetic disorder can be used in 
simulations of structural phase transitions in a vicinity 
of magnetic Tc"^*^''^' . Still, to the best of our knowledge 
all the applications of the DLM model so far neglected 
the effect of lattice vibrations. 

On the other hand, the importance of lattice dynamics 
for an accurate description of thermodynamic properties 
of materials is well recognized by now ' . In particu- 
lar, its consideration can be essential for a treatment of 
lattice stabilities'""^', heat capacities'''' ', and equations 
of state''' of solids. The importance of lattice vibrations 
should be fully recognized in the paramagnetic state of 
magnetic materials, as it occurs only at elevated tempera- 
tures. However, simultaneous treatment of the magnetic 
disorder, inherent to the paramagnetic state, and lattice 
vibrations represent a truly challenging task. 

State-of-the art treatments of lattice vibrations are 
based either on (quasi-) harmonic calculations of the 
phonon dispersion relations or on molecular dynamics 
simulations'' ' . For magnetically ordered materials these 
techniques can be applied straightforwardly. However, in 
the presence of magnetic excitations this will not work. 
In particular, in the paramagnetic state at high temper- 
atures the relevant magnetic excitations are associated 
with spin-flips. Their characteristic time scale can be es- 
timated by the spin decoherence time tdc- Spin-dynamics 
simulations of the spin autocorrelation function in bcc Fe 
above T^'' show that is of the order of 20-50 fs. In 
materials with lower Tc it should be larger by approxi- 
mately a corresponding factor, because both Tc and the 
velocity of the propagation of the local moments are re- 
lated to the strength of the exchange interactions. At 
the same time a typical MD run should be carried out 
for at least 3-5 ps, as dictated by the inverse of the De- 
bye frequency (~ 10"'^^ s). This means that magnetic 
configurations should change often during the MD run. 
Simultaneously, a typical MD time step is of order of 1 fs, 
which is still much smaller than tdc- Thus, the magnetic 
degree of freedom is slow on the time scale relevant for 
the determination of temporal evolution of a particular 
atomic configuration, but fast on the time scale relevant 



for a proper exploration of the phase space of atomic con- 
figurations. Therefore, the adiabatic decoupling between 
magnetic and vibrational degrees of freedom cannot be 
applied, and they should be treated within one single 
framework. Similar arguments can be used to question 
a validity of lattice dynamics studies for paramagnetic 
materials based on quasi-harmonic approximation. Per- 
haps the most consistent approach to the analysis of spin- 
lattice interactions at finite temperature would be to ap- 
ply a combination of molecular dynamics with ab ini- 
tio spin dynamics'^" or with DMFT. However, at present 
such calculations are hardly feasible in practice. 

Within our approach we describe the paramagnetic 
state of a system within the disordered local moment 
picture. In this approach, local moments exist at each 
magnetic site of a system (in our case, at Cr sites in 
CrN) and are commonly thought to fluctuate fairly inde- 
pendently. Thoughtful discussions of the DLM model can 
be found in Ref.s. '.28,30-32^ r^j^^ ^^^^^^ ^j^^ j^j^^^ 

proach in the many-body lattice models like the Hubbard 
or s-f exchange "Kondo lattice" model is discussed in 
Ref. , where it is argued that though in a complete the- 
ory the charge and spin fields are dynamically fluctuating 
both in space and time, a "static" DLM approximation, 
where one neglects the dynamics of the fluctuations cap- 
tures an important part of the correlations. In the DLM 
a correlated system is described in terms of a pseudo- 
alloy of spin up and spin down components. Combined 
with the coherent potential approximation (CPA) it be- 
comes equivalent to the "Hubbard HI" approximation'"' 
for the original many-body problem, and it is used with 
success in many applications. 

However, the CPA is applicable for the description of 
a substitutionally disordered system with atoms at the 
sites of an ideal underlying crystal lattice , and there- 
fore cannot be used for treatment of lattice dynamics at 
finite temperatures. In a previous work , we took one 
step towards the simultaneous modeling of magnetic and 
vibrational finite temperature effects by suggesting two 
alternative supercell implementations of the DLM calcu- 
lations. In the first, a specific collinear distributions of up 
and down magnetic moments arranged to minimize the 
spin correlation functions were used, in line with the spe- 
cial quasirandom structure (SQS) methodology. In the 
second, a magnetic sampling method (MSM) was pro- 
posed. In the MSM, the energies of a number of ran- 
domly generated magnetic distributions were calculated 
and their running average was taken as the potential en- 
ergy of the paramagnetic sample. In Ref. it was shown 
that for CrN MSM calculations are converged already 
for 40 different magnetic distributions, and the two ap- 
proaches, the SQS and MSM give almost identical results. 

The SQS approach makes use of the fact that in a static 
picture with all atoms fixed on ideal lattice points, the 
description of a spatial disorder between the local mo- 
ment orientations is a good approximation to model the 
energetics of the combined space and time fluctuations of 
magnetic moments in a real paramagnet. Unfortunately, 



3 



if the vibrations of atoms are to be included, one needs to 
go beyond the fixed magnetic state described by the SQS. 
The reason is that if a magnetic state is fixed in time one 
would see artificial static displacements of atoms off their 
lattice sites due to forces between the atoms with differ- 
ent orientations of their local moments and with different 
local magnetic environments. In the CrN case those are 
likely to be quite large due to the magnetic stress dis- 
cussed in Ref. '. In a real paramagnet, due to the time 
fluctuations of the local moments, this effects should be 
at least partially averaged out and suppressed depending 
on the time scales of the spin fluctuations and atomic 
motions. 

The MSM could in principle be used to obtain the 
adiabatic approximation where the magnetic fluctuations 
are considered to be instantaneous on the time scales of 
atomic motions. This approximation would be obtained 
if the forces acting on each atom were averaged over a 
sufficient number of different magnetic samples during 
each time step of a molecular dynamics simulation. The 
obvious drawback in this approach is that a large num- 
ber of calculations needs to be run in parallel leading to 
an increase, a factor 40 in our case, in computational 
demands. Furthermore, as stated above it is not at all 
clear that this adiabatic approximation is motivated in 
any system. However, the MSM gives us a very good 
starting point for the implementation of the DLM pic- 
ture in a MD framework. 



B. Disordered local moments molecular dynamics 

In this work we introduce a method for molecular dy- 
namics simulations of paramagnetic materials within the 
traditional ab-initio MD framework. Starting from the 
DLM idea of a spatial disorder of local moments, we also 
change the magnetic state periodically and in a stochastic 
manner during our MD simulation. In this way we deal 
with a magnetic state that does not show order either on 
the length scales of our supercell, or time scale of our sim- 
ulation. We make an approximation that the magnetic 
state of the system is completely randomly rearranged 
with a time step given by a spin flip time {Atgf), and 
with a constraint that the net magnetization of the sys- 
tem should be zero. Hence to simulate a paramagnetic 
system with a spin flip time At^f, we initialize our cal- 
culations by setting up a supercell where coUinear local 
moments are randomly oriented and the total moment of 
the supercell is zero, and run coUinear spin-polarized MD 
for the number of MD time steps (AImd) corresponding 
to the spin flip time, that is for Atsf / AImd time steps. 
Thereafter the spin state is randomized again, while the 
lattice positions and velocities arc unchanged, and the 
simulation run continues. 

Here it is worth to point out that besides the treatment 
of the many-body effects important for the description of 
the paramagnetic state at the DLM-LSDA level, or as will 
be discussed below for the present case: DLM-LSDA-I-U, 



we introduce several additional approximations. In par- 
ticular, we neglect effects due to non-coUinear orienta- 
tions of the local magnetic moments. This, however, is 
justified for the paramagnetic state well above the mag- 
netic transition temperature' . Note also that magni- 
tudes of the local magnetic moments are allowed to vary 
as dictated by the self-consistent solution of the elec- 
tronic structure problem at each step of MD simulation. 
At the same time, we substitute the true spin dynamics 
with instantaneous modification of the sample magnetic 
structure with time steps Atgf. Here we follow Ref.^'^ 
and make use of the physical picture that the simulated 
system, although ergodic, does not cover its phase space 
uniformly in time. In the DLM model one assumes that it 
gets stuck for long times, of the order of At^f, near points 
characterized by a finite moment at every site pointing in 
more or less random directions and then moves rapidly 
(in our case instantly) to another similar point. The 
states of temporarily broken ergodicity are characterized 
by classical unit vectors, ej, assigned to each site i and 
giving the direction of the magnetization averaged over 
the spatial extent of the i-th site in the supercell and the 
time Atgf . The motion of temporarily broken ergodicity 
is mainly characterized by changes in the orientational 
configuration of the moments. 

Note that in Ref. ' ' the magnetic degree of freedom was 
related to an inverse spin- wave frequency tgu; ~ 1 /ujsw 
100 fs, which represents the dominating magnetic excita- 
tion at low temperatures. However, in the paramagnetic 
state at high temperatures the relevant magnetic exci- 
tations are associated with spin-flips rather than with 
spin waves. Thus the relevant time scale is better char- 
acterized by the spin decoherence time At^c rather than 
by the inverse spin-wave frequency. As we pointed out 
above, the latter was estimated to be of the order 20-50 fs 
in bcc Fe above Tq ■ For CrN, with a Tjv around room 
temperature and probably with weaker exchange inter- 
actions, we expect that tdc could be somewhat larger. 

However, our procedure makes it possible to model a 
paramagnetic system for any particular time scale of the 
spin dynamics. In fact one can span the whole range be- 
tween the two adiabatic approximations: from the frozen 
magnetic structure to magnetic configurations that rear- 
range instantaneously on the time scales of each atomic 
motion during the MD run. Of course, the appropriate 
value of this parameter needs to be found with real spin 
dynamics calculations or taken from experiments. In this 
paper, we study a range of different spin flip times and 
their consequences for the obtained structural and ther- 
modynamic properties of CrN. 



C. Calculational details 

All our flrst-principles calculations in this work is per- 
formed using the projector augmented wave method as 
implemented in the Vienna Ab-initio Simulation Pack- 
age (VASP) The electronic exchange-correlation 
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effects are modeled using a combination of the Local 
Density Approximation with a Hubbard Coloumb term 
(LDA+U) using the double-counting correction scheme 
suggested by Dudarev et al:'''. The value of the ef- 
fective U {U'^^f ^ U - J) apphed only to the Cr 3d- 
orbitals is taken as 3 eV, found to be suitable in com- 
parison with several experimentally measured structural 
and electronic properties in Ref. . 

Our simulation box, both in the simulation of the cubic 
and orthorhombic phases, contains 32 Cr and 32 N atoms 
arranged in a supercell of 2x2x2 conventional cubic unit 
cells. In the orthorhombic case, the primitive vectors of 
the supercell are tilted and scaled in line with the re- 
sults of a structural optimization of this low temperature 
antiferromagnetic phase. 

The plane wave energy cut-off was set to 400 eV. We 
used a Monkhorst-Pack scheme " for sampling of the 
Brillouin zone using a grid of 2x2x2 k-points. To check 
the accuracy of the potential energies and pressures a 
selection of configurations are chosen out of the MD sim- 
ulation run and recalculated with a higher accuracy. The 
error arising from for the k-point sampling is relatively 
constant with a shift of about 35 meV and a standard de- 
viation of less then 2 meV. Hence the relative potential 
energies that are calculated have a high accuracy. The 
pressures also have a small constant shift of about 0.2 
GPa with a standard deviation less then 0.1 GPa when 
the electronic structure calculations are converged with 
respect to the k-point mesh. 

In order to control the temperature of the simulation, 
avoid artificial energy drift and to minimize the influence 
of the particular choice of initial magnetic and lattice con- 
figurations in the simulations we use a Nose thermostat 
The values of the bulk modulus, K^, have been deter- 
mined by fitting our calculated pressure and volume data 
to the Birch-Murnaghan equation of state"'''"''". 

P = 3Xo/b(1 + VeT'^ (l + 2/3(< - 4)/^^) (1) 

where is the bulk modulus and Kq is the derivative 
of Kq with respect to pressure. The fulerian strain, is 
defined as fs = l/2[(Vo/V")2/3 - 1], where V and Vq are 
the volume respectively equilibrium volume. 

III. APPLICATION TO CrN 

A. The potential energy 

From the MD calculation we extract the potential en- 
ergies of CrN. As can be seen in Fig. 1, the potential 
energy of the system is well conserved. To investigate 
the infiuence of the spin fiip time, the potential energy of 
CrN is calculated for several Ais/. The results are shown 
in Fig. 1. 

In Fig. 2 these potential energies are collected and 
shown relative to the potential energy of the calculations 



with shortest Atg/, 5 fs. There is a clear shift in poten- 
tial energy of about 10 meV from the simulations with 
the shortest spin flip times of 5 fs to the longest of 100 
fs. This can be compared with the total energy reduction 
due to static relaxations of 15 meV that we get by using 
the SQS approach treating the magnetic state as frozen 
in time. We note that for the lower values of At^j, corre- 
sponding to fast spin decoherence, there is a plateau were 
the potential energy is only weakly influenced by At^/. 
However, between spin flip times of 15 to 50 fs, there is 
a considerable change in potential energy. Of course, the 
energy scale should be material specific. We suggest, as 
a quick test of the importance to consider this effect, a 
calculation of relaxation energies of a paramagnetic sys- 
tem using the SQS approach'' with a fixed magnetic state 
through the relaxation. The obtained relaxation energy 
should correspond to an upper limit on the potential en- 
ergy dependence on Ai^/. 
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FIG. 1. Potential energy of cubic paramagnetic CrN as a 
function of simulation time calculated at 300 K using DLM- 
MD method. Shown are the results obtained with a spin flip 
time of 10 fs and 100 fs, as well as with a static magnetic 
state. Results for conventional AIMD simulations carried out 
for CrN in the orthorhombic antiferromagnetic ground state 
are also shown for comparison. The potential energy is stable 
and well converged as can be seen by the included running 
averages. 



B. Pair distances 

In order to analyze the difference between the pro- 
posed DLM-MD simulations and magneto-static MD in 
more details, an investigation of the local environment 
of the different atoms is carried out, especially the Cr 
- Cr metal nearest neighbor distances. In Fig. 3 his- 
tograms are shown of all the Cr - Cr nearest neighbor 
distances. These are also separated into ft; ii and |t 
pairs. Hence we can see the effect of the magnetic state 
on the distribution of pair distances. In Fig. 3a the Ats/ 
is very short, 10 fs, hence the atoms do not have time 
to adjust their positions for the current orientation of lo- 
cal magnetic moments and we do not see any difference 
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FIG. 2. Potential energy shift calculated for paramagnetic 
CrN as a function of the spin flip time Ataf. The shortest 
spin flip time of 5 fs is taken as reference. 
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FIG. 3. Histogram of the Cr - Cr nearest neighbor distances 
for Cr atoms with parallel (solid line) and antiparallel (dashed 
line) orientations of local magnetic moments obtained from 
DLM-MD simulation for the cubic paramagnetic phase at 300 
K. a) Afs/ =10 fs. b) Afs/ =100 fs. seen, c) Static magnetic 
state, d) The orthorhombic antiferromagnetic phase of CrN 
calculated with conventional AIMD for comparison. 



comparison. Here the ft: 4-4- a-nd ti: it pairs of magnetic 
moments are arranged in an ordered way, see e.g. Fig. 
4 in Ref.'', that allows for maximal relaxation of atomic 
coordinates in combination with a structural relaxation 
of the unit cell, giving rise to a large separation between 
the two different kinds of pairs. 

A possibility of statistical correlations between the 
atomic distances and the orientation of atomic moments 
also in a dynamically changing paramagnetic phase is in- 
deed an intriguing thought experiment. Although we can 
not rule out its existence from principal considerations, 
to the best of our knowledge it has never been reported 
in experiments. However, we note that our two main 
approximations in the present DLM-MD, the usage of 
coUinear moments and the temporarily broken ergodicity 
of the DLM approach, are likely to introduce inaccura- 
cies that exaggerate those local spin-lattice correlations 
when a slow spin dynamics is modelled. Therefore, when 
the here suggested method is used, a smaller value of 
A^s/ , corresponding to the absence of differences in dis- 
tances between atoms with parallel and antiparallel local 
moments, like in Fig. 3a, should be recommended. The 
presence of an energy plateau in Fig. 2 seems to indicate 
that this should be a reasonable approach. 



C. Bulk modulus of CrN 




Pressure (GPa) 



in distances between the ttiii and the tiiit pairs. In 
Fig. 3b, the spin flip time is increased to 100 fs and now 
the atoms have had sufficient time to move towards the 
energetically preferential positions. Consequently, a shift 
in pair distances between ftj and the fi, it is evident. 
Fig. 3c is obtained with the same orientation of local 
moments during the whole MD run. Here we also see a 
splitting in the pair distances between the tti ii and the 
ti, it pairs, which is of the same order as in the previ- 
ous case. Hence, 100 fs between the re-arrangement of 
magnetic configurations is long enough for the atomic nu- 
clei to adjust considerably their positions in the supercell 
to the given magnetic configuration. In the last figure. 
Fig. 3d, the pair distances are shown for the low temper- 
ature antiferromagnetic orthorhombic ground state for 



FIG. 4. Volume as a function of pressure for the cubic and 
orthorhombic phase from MD simulations at 300 K. The equa- 
tion of state for the cubic phase is calculated using a spin flip 
time of 10 fs. The calculated volumes are normalized with 
the calculated equilibrium volume of the cubic phase, and the 
experimental points' with the measured equilibrium volume 
of the cubic phase. The inset shows the dependence on spin 
flip time for the calculated equation of state 

Our main goal with this work is to study the equation 
of states, and in particular the bulk modulus of paramag- 
netic CrN which has recently been discussed in the litera- 
ture' Using our DLM-MD approach we are able to cal- 
culate volume as a function of temperature and pressure 
for both the paramagnetic cubic and the antiferromag- 
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TABLE I. Calculated bulk modulus, Ko, of CrN in the or- 
thorhombic antiferromagnetic and the cubic paramagnetic 
phases obtained at ambient pressure and 0, 300, and 1000 
K respectively. 
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netic orthorhombic phases. In the former case we also 
investigate if there is an impact of the value of the spin 
flip time parameter on the equation of state. Thus we 
are able to investigate how both the dynamical change of 
magnetic configurations in the paramagnetic state and 
the lattice vibrations, neglected in previous theoretical 
works but of course present in the experiments, impact 
on the compressibility. Figure 4 shows the calculated 
volume as a function of pressure for the two phases at 
300 K and compare them to the experimental measure- 
ments by Rivadulla et al. ' . One sees very good agree- 
ment between theoretical and experimental equations of 
state. In particular the relative shift in volumes between 
the two phases is reproduced within the measured er- 
ror bars. The calculated slope of the orthorhombic phase 
agrees well with the measured values for this phase where 
the measurement is done over a large pressure range. The 
inset in Fig. 4 shows the influence of the spin flip time on 
the volume versus pressure curves in paramagnetic cubic 
CrN. A change in Atgj introduces a small shift of the 
volumes, but does not influence the slope of the curves. 

Our results confirm the possibility of a pressure in- 
duced phase transition from the cubic paramagnetic to 
the orthorhombic antiferromagnetic phase due to the 
slightly smaller volume of the latter, in line with previous 
investigations. Importantly, as can be seen in Table I and 
from the slopes of the curves in Fig. 4, the bulk modu- 
lus is found to be very similar between the two phases. 
This is the case both at 300 K, 1000 K, and in the static 
K calculations. The calculations of the orthorhombic 
low temperature phase at 1000 K is of course not of rel- 
evance for any comparison with experiments, but is in- 
cluded to show with certainty that temperature induced 
vibrations is not influencing the difference in bulk mod- 
ulus between the phases. At T=300K and P=OGPa we 
find = 290 GPa while if^^*^ = 286 GPa. This 

gives an insignificant difference of 4 GPa, far from the 
collapse of 25% or 85 GPa suggested in Ref. to follow 



the transition from cubic to orthorhombic structures. A 
variation of the time between the rearrangement of the 
magnetic configurations do not influence the value of the 
bulk modulus in any appreciable way. Thus, explicit con- 
siderations of temperature induced magnetic fluctuations 
and lattice vibrations do not change the main conclusions 
from previous works ' : There is no theoretical support 
for a collapse of the bulk modulus of CrN upon the pres- 
sure induced phase transition. 

IV. SUMMARY 

We present a method for calculation of thermodynamic 
properties of magnetic materials in their high temper- 
ature paramagnetic state. We use ab-initio molecular 
dynamics and simulate the paramagnetic state with dis- 
ordered non- vanishing local magnetic moments. Random 
configurations of the local moments in the simulation cell 
are switched at predetermined time intervals. Hence we 
can capture the influence of the dynamically disordered 
magnetic state on the lattice dynamics as it develops 
during the simulation. We apply this method to CrN 
which is known to have a strong interaction between the 
magnetic state and the lattice. We find that there is a 
connection between how fast the local moments are al- 
lowed to flip and the calculated potential energy. If the 
spin flip time is short, ^ 10 fs, the lattice do not have 
time to respond, but if the spin flip time is increased 
to about 100 fs then the atomic positions start to show 
clear relaxation effects. We apply this disordered local 
moments molecular dynamics method to the calculation 
of the equation of state of paramagnetic cubic CrN and 
compare with calculations for the orthorhombic antifer- 
romagnetic phase and with experiments. In particular we 
calculate the debated bulk modulus and find that there 
is only a very small difference, and definitely no collapse, 
in Kq between the orthorhombic antiferromagnetic phase 
and the cubic paramagnetic phase. 
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